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1 Introduction 

Fractional calculus has found motivations in a growing area concerning gen- 
eral stochastic phenomena. These include the appearance of alternative dif- 
fusion mechanisms other than Brownian, as well as classical and quantum 
mechanics formalisms including dissipative forces, and therefore allowing an 
extension of the quantization schemes for non-conservative systems [I]- 11]. In 
particular it is interesting to study the time fractional Schrodinger equation, 
which describes the time irreversibility present in all processes in nature. 
This inhomogeneity in time is particularly important at the quantum level, 
where damping is inherent to all non equilibrium evolution. Our aim here is 
to extend a previous study about this equation, introduced in [3J. With this 
purpose we use an analytical definition of fractional derivative [5], where we 
show that it is possible to obtain a general solution for the time fractional 
equation, for any complex value of the derivative index. Furthermore the 
associated Green functions can be evaluated in a straightforward way. This 
analytical generalization is possible within the framework of Ultradistribu- 
tions. 

Ultradistribution's theory (ref.[TTJ [12]) proofs to be a mathematical ad- 
equate theory to approach several branches of theoretical physics, helping 
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to overcome certain common subtleties arising from the formal handling of 
analytical formulas. For instance, using this tool in ref. [7] we have realized 
the mathematical treatment of higher order field theories, and in ref. [8] we 
have solved the problem of normalization of Gamow states or resonances in 
quantum mechanics. 

Ultradistributions have the advantage of being representable by means of 
analytic functions. So that, in general, they are easier to work with them. 

They have interesting properties. One of those properties is that Schwartz 
tempered distributions are canonical and continuously injected into tem- 
pered ultradistributions and as a consequence the Rigged Hilbert Space with 
tempered distributions is canonical and continuously included in the Rigged 
Hilbert Space with tempered ultradistributions. 

Therefore the properties of ultradistributions are also well adapted for 
their use in fractional calculus. In this respect we have shown that it is 
possible (ref. [5]) to define a general fractional calculus with the use of them. 

This paper is organized as follow: 
In section 2 we define the fractional Schrodinger equation for all v complex 
with the use of the fractional derivative defined via the theory of tempered 
ultradistributions. In section 3 we solve this equation for the free particle and 



4 



give three examples:^ = 1 /2, v = 1 and v = 2. In section 4 we realize the 
treatment of the potential well and we analyze the cases y = 1 /2, y = 1 and 
"v = 2. In section 5 we study the Green fractional functions for the free parti- 
cle in three cases: the retarded Green function, the advanced Green function 
and the Wheeler-Green function. As an example we prove that for y = 1 
these functions coincide with the Green functions of usual Quantum Me- 
chanics. In section 6 we discuss the results obtained in the previous sections. 
Finally we have included three appendixes: a first appendix on distributions 
of exponential type, a second appendix on tempered ultradistributions and 
a third appendix on fractional calculus using ultradistributions. 

2 The Fractional Schrodinger Equation 

Our starting point in the study of the fractional Schrodinger equation is the 
current known Schrodinger equation: 



2m 



a^(t,x) + v(x)i|)(t,x) 



(2.1) 



According to ref. [3], (12. ip can be written as: 



iT p 9 t t|j(t,x) 



V(x) 



■i|>(t,x) 



(2.2) 



2m 
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where L P = a/ Gh/c 3 , T p = a/ Gh/c 5 , M p = ^Jhc/G and E p = M p c 2 , and G 
is the gravitational constant. 

If we define N m = m/M p and N v = V/E p we obtain for fT2T2|) 



Li 

iT p 3 t i|)(t,x) =-^9^(t,x)+N v il)(t,x) (2.3) 



< in 



By analogy with ref.p] we define the fractional Schrodinger equation for all 
v complex as: 

Li , 

(iT p ra^(t,x) =-^-a^(t,x) + N v ^(t,x) (2.4) 

where the temporal fractionary derivative is defined following ref.[5] (see 
Appendix III) 



3 The Free Particle 

From (12.41) for the free particle the fractionary equation is: 

(ia t )>(t,x) + -i-a^( t ,x) = o (3.i) 

By the use of the Fourier transform (complex in the temporal variable and 
real as usual in the spatial variable) the corresponding equation is (see Ap- 
pendix II and ref.[5]) 

- T^k 2 I <p(ko, k) = b(ko, k) (3.2) 



2T p -N m 
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whose solution is: 



b(ko,k) 



"0 2T^N, 



-k 2 



(3.3) 



and in the configuration space (anti-transforming) 



ib(t,x) = () 



a(ko,k) 



r -oo 2T^N, 



-k 2 



. e -i(k t + kx) dk dk 



(3.4) 



where: 



a(ko,k) 



b(k ,k) 
An 2 



As a(ko, k) is analytic entire in the variable ko we have for it the development: 



a(ko,k) = Y_ am(k)k» 

m=0 

Inserting this development in (13 .41) we obtain 



iKt,x) = Y_° 



m=0 



a m (k)k, 



L e -i(k t+kx) dk dk 



-k 2 



(3.5) 



r -oo 2T^N r 

In ref.j6] is given the Laplace transform of the generalized Mittag-Leffier's 
function 



e-^- ] ^(x a z) dx 



,a-|3 



|J. a — Z 



(3.6) 



and as a consequence 



v— m— 1 1 



-v,v— m 



[(x + iO) v e- l T-z]e lkx dx = 



k v -z 



(3.7) 



with 3(k) > 0, and 



e -if (v-m) (x + to)^- m -'E V)V _ m [(x + i0) v e-^z]e l,cx dx 



* v— m— 1 



ikx 



k v -z 

with 3(k) < From which we obtain for the complex Fourier transform 



(3.8) 



Hp(k)] 



e -i f iv-mj( x + to)^™-- 1 ^ v _ m [(x + i0) v e-^z]e to dx 



ikx 



H[-Df(k)] 



e -x Tl v-m]^ x + iO) v - m -'Ev, v _ m [(x + i0) v e^-z]e lkx dx 



k v -z 

and then we have for the inverse complex Fourier transform: 



1 

In 



O: 



k v -Z 



e -t- ( v-m, (x + to)^-™--iE V)V _ m [(x + i0) v e-^z] (3.10) 



_ e - ikx dk = 



Thus we have from ( 13. 5p 



x|j(t,x) =^e- i 5 (v - m) (t + i0) 



v— m— 1 



X 



m=0 



a m (k)F v , v _ m [(t + i0) v e-^w]e- ikx dk 



where we have made the rescaling ^a m (k) — > a m (k) and 



2 1,2 



w = 



14k 



2T V N 



(3.9) 



(3.11) 



When a(k , k) = a(k) formula (13.111) simplifies: 



t|)(t,x) = e-^ iy) {t + iO) 



\y — i 



a(k)E ViV [(t + iO) y e- iI! rw]e- ikx dk (3.12) 



When we consider the solution only for t > 0, i.e. writing il>(t, x)H(t) in 
(13. ip . we obtain the result of ref.[3]. In fact, the equation for this case is: 



(i3t) v -'{[(i3t)M>(t,x)]H(t)}- 



l " a^(t,x)H(t)] = o 



2N J y 



(3.13) 



With the initial condition ib(0, x) = ih(x) we have for (I3.12p 



\1 . 

v 8j[il)(t,x)H(t)] = (3.14) 



2N T v 



Let ij)(ko,k) and ih(k) the Fourier transforms of ib(t, x)H(t) and x|)(x) re- 
spectively. Then ( 13.141) transforms into: 

kJ<j>(lco,k)-ilcJ- 1 $(k)HDf(1co)]-w$(ko,lc) = (3.15) 



whose solution is 

$(k ,k) = i 
Using again ( 13. 6 h we have: 



-S_L{j5(ic)H[a(ko)] 

k^-w 



(3.16) 



Hp(k)] 

and anti-transforming 



E v [(-ix) v z]e lkx dx = H[3(k)]i 



k v -z 



(3.17) 



1 

In 



.v-1 



d>H[3(k)]i- e~ ikx dk = H(x)E v [(-ix) v z] 

k v — z 



(3.18) 



Then taking the inverse Fourier transform of (I3.16P we obtain 



H(tMt,x) = 



H(t) 
In 



Ej(-it) v w]iL.(k) dk 



(3.19) 



Thus, for t > we have: 



i|)(t,x) 



1 

2tt 



E v [(-it) v w]iL)(k) dk 



(3.20) 



which is the result obtained in ref.[5]. 

We proceed to analyze the solutions of ( 13 .4p for some typical cases of the 
derivative index v in the following section. 



Examples for free particle fractionary time evolution 



As a first example we consider the case v = 1 /2 



Let a be given by: 



a 



4 



2T p 'N m 



From ( 13.4ft we obtain 

i|)(t,x) 

or equivalently: 

i|)(t,x) 



a(k , k) _ t 



e -i ( k t+kx) dk dk 



r -oo k o ~ al<2 



oo oo 

-i(a 2 k 4 t+kx) 



a(k)e" 



J dk- 



a(ko,k) 



(3.21) 



(3.22) 
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1 



1 



(k + iO) 2 - <xk 2 (ko - iO) z - <xk 2 



-i(k t + kx) dko dk 



where: 



(3.23) 



a(k) = -4mcck 2 a(cc 2 k 4 ,k) 



With some of algebraic calculus we obtain for OIL 



a(k)e 



-i(tu 2 t+kx) d i < ._ 



^l^e^^-^dko dk 
k + cu 2 



-oo 



with: 



tu = ak 2 



(3.24) 



and where we have made the re-scaling: 

i 

-2ik 2 a(-k ,k) — > a(ko,k) 

The first term in (13 .8p represent free particle on-shell propagation and the 
second term describes the contribution of off-shell modes. 

As a second example we consider the case v = 1 . 
In this case ( 13 .4[) takes the form: 



4)(t,x) = o 



k — w 



(3.25) 



r — oo 
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Evaluating the integral in the variable ko we have: 



l|)(t,x) 



a(k)e- i(u,t+kx) dk 



(3.26) 



where a(k) = — 2ma(u), k). Thus we recover the usual expression for the 
free-particle wave function. 



Finally we consider the case y = 2. For it we have 



\b(t,x) = o 



a(ko,k) ^_ i(kot . 



-kx) 



r — oo 



kg-CU 2 



dk dk 



(3.27) 



After to perform the integral in the variable ko we obtain from (13. 27ft : 



ib(t,x) 



a{k)e- i[wt+kx] + b+(k)e i(a,t+kx) dk 



(3.28) 



with a(k) = — 2ma(cu,k) and b + (k) = — 2ma(— cu, — k) 

4 The Potential Well 

We consider in this section the potential well. The fractionary equation for 
a particle confined to move within interval < x < a is: 



(iT P raM>(t,x) 



li , 
P 3^(t,x) 



2N 



(4.1) 



To solve this equation we use the method of separation of variables. Thus if 



we write: 



xMt,x) = ih (t)ib 2 (x) 



(4.2) 
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As is usual we obtain: 



i|>i(t) 

Then we conclude that "4>2(x) satisfies: 



= A (4.3) 



3^ 2 W + ^ib 2 (x)=0 (4.4) 
p 



The solution of (14. 4 h is the habitual one: 



V a 

with: 



ib 2n (x) = b n sin ( — x ) (4.5) 



An = ^(^ 2 (4.6) 



2N 



m 



and the boundary conditions satisfied by 4>2n( x ) are: 

As a consequence of ( 14.3ft .( 14751) and (14.61) the Fourier transform ibi(ko) of 
ibi (t) should be satisfy: 

(l^-w n )$ ln (ko) =0 ( 4 -7) 
(w n = A n /Tp ) whose solution is: 

^Inl^o) = 7^ 4.8 
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Therefore the final general solution for \Jj> (t, x) is: 



ib(t, x) = ^ sin (— — x) o 



n=l 



V a 



dko 



(4.9) 



where we have defined: 



a n (k 



bnCuCko 
In 



which is an entire analytic function of ko. As a consequence we can develop 
a n (ko) in a power series already of k = 



a n (k ) = ^ a nm k, 



and obtain for i|)(t, x) the expression: 



i|)(t,x)= Y" a nm sin(— x)o - e-^dko (4.10) 
. . V a / . kg — w n 



n=1 ,m=0 

According to (I3.10p we have then: 



i|)(t,x)= ^ a 

n=1 ,m=0 



n7i 



nm sm( — x)e " w (t + iOP ' x 



\v— m— 1 



V a 



If we select a n (k ) = a n independent of the variable ko, we obtain: 



00 

i|>(t,x) = Y_ Q n sin (-- x ) e _i 5 (v) (t + i0) 



\y — l 



n=l 



V a 



x 



(4.11) 



(4.12) 
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When we consider the solution only for t > 0, i. e. ib(t, x)H(t), we recover 
the result given in ref.[3], since the equation for this case is: 

(i9 t r 1 {[(i3 t )^(t,x)]H(t)}+ -il— 85[iKt,x)H(t)] = (4-13) 

With the initial condition i]>(0, x) = ih(x) we have for (14. 13j) 

(ia0^(ia t )[iJ)(t,x)H(t)]-iti)(x)5(t)}+-i- 9^(t,x)H(t)] = (4.14) 

Purposing again 

ip(t,x) = xbi (t)\b 2 (x) 

we obtain two equations: 

j;(id t r-\(id t ) [ip, (t)H(t)] - (0)5(t)} - A^! (t)H(t) = (4.15) 

a^ 2 (x) + ^x|; 2 (x)=0 (4.16) 
p 



Solution of (I4.16P is again: 



V a 

with: 



ib 2n (x) = bn sin ( — x ) (4.17) 



*» = ™-(VJ (418) 



2N 



m 



Using the Fourier transform for (I4.15P we obtain: 

T;k^ (k ) - ikr>! (O)HCDf(ko)] - V$1 (ko) (4.19) 
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whose solution is: 



p 



Icq - w n 



and then 



tl) 1 (t)H(t) = ^^ct)H[j(k )l 
27tT p v 



k^ 1 
-e- ikot dk 



-w„ 



•0 



In the configuration space (I4.2ip reads: 



^(t)H(t) = ^lH(t)E v [(-it) Y w n ] 



Then we have for ip(t,x): 

oo 

xp(t,x) = ^ a n sin (-^-*) ^[(-it^wj 



n=0 



(4.20) 



(4.21) 



(4.22) 



(4.23) 



Examples for Potential Well fractionary time evolution 



As a first example we consider the case "V = 1/2. For it the solution (14. 9 p 
takes the form: 



n=1 



ib(t, x) = sin ^ — x 



rajlcoje-^ 

° — I 

•J k 2 - w n 



dk 



(4.24) 



or equivalently: 



oo 

xp(t,x) = ^ a n sin (— x) e" iw ^+ 

n=1 



n=1 



/ TL7T 
sm X 



V a 



(k + iO) 2 — w n (ko — i0)2 — w r 



x 
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a n (k )e- ikot dk 



(4.25) 



After performing some algebra we have for ( 14.111) the expression: 

oo 

ib(t,x) = Y_ Q nsin (™*) e" iw ^+ 



n=l 



/ TL7T 

sin y — x 



a n(k)) p -i1c t 



dk 



(4.26) 



k + 

Analogously as before, the second term in f)4.12p represents of-shell stationary 
modes. 



As a second example we consider v = 1 . In this case: 

if. \ v 1 • fnn \ f Q n (ko)e- ikot 

ip[t, xj = > sin — x () — dk 

\ a J ] k - w n 

TL— I p 

Performing the integral in the variable ko we have: 



oo 

lL>(t,x) = Y_ a ^ sin (~ X ) C ^ 
n=1 



Which is the familiar general solution for the infinite well. 



Finally for v = 2: 



ib(t, x) = sin ^-^x 



o 



n=l 



[ Mkoje-^ 



k^ -w n 



dko 



and after to compute the integral: 

oo 

xp(t, x) = Y_ s in (— *) (due- 1 ^ + b n e +i ^ t 
with a n = a n {y/w^) and b+ = a n ( 



(4.27) 



(4.28) 



(4.29) 



(4.30) 
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5 The Green Function for The Free Particle 

As other application that shows the generality of the fractional calculus 
defined with the use of ultradistributions, we give the evaluation of the Green 
function corresponding to the free particle. Let (3 be defined as: 

p2 = rF!r <"> 

Then G(t — t',x — x ) should be satisfy the equation: 

(i9 t ) v G(t-t',x-x') + |3 2 3 2 G(t-t',x-x') = 5(t - t')6(x - x') (5.2) 

As G is function of (t — t',x — x') it is sufficient to consider G as function of 
(t,x): 

(ia t ) v G(t,x) + |3 2 9 2 G(t,x) = 5(t)5(x) (5.3) 
For the Fourier transform 6 of G we have: 

flcj - |3 2 k 2 )G(k ,k) = Sgn[ ^ ko)] + a(k ,k) (5.4) 

where a(ko, k) is as usual a rapidly decreasing analytic entire function of the 
variable k . Selecting: 

a(k ,k) = - 

we obtain the equation for the retarded Green function: 

(k£ - (3 2 k 2 )G Tet (k ,k) = Hp(lco)] (5.5) 



and then: 



G ret (t,x) = ^\> 



r — oo 



If we take: 



a(k ,k) 



we obtain the advanced Green function: 



Gadv(t> xj 



An 2 



\ H[-J(Tcq)] i(kntH 

k^ - |3 2 k 2 

r — oo 



kx) dk dk 



For the Wheeler Green function (half advanced plus half retarded): 



G w (t,x) = - [G Qd v(t,x) + G ret (t,x)] 



we have: 



Gw(t,x) = i ^ 



r — oo 



f Sgnp(ko)] ^ 
k y - |3 2 k 2 



(k t+kx) 



dk dk 



Using ( 13. 7p . ( 13. 8 p with m = we obtain: 

oo 

1 



G Te t(t,x) = ^^t;- 1 



G Q dv(t, x) 



271 



e -i7t(|-1) t Y-1 



G w (t,x) = -Le- iZ ¥ i (t + iOr- 1 

471 



t v ,v[e^-(t + i0) v (3 2 k 2 ]e- tkx dk 
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Example for the free particle Green function 



When we select v = 1 we obtain the usual Green functions of Quantum 



Mechanics. For example for G Tet we have: 



G Tet (t,x) = 4^9 



f H[IJ(k, 



k - |3 2 k 2 



OJJ -i(k t+kx) 



dk dk 



(5.13) 



r — oo 



or equivalently: 



G ret (t,x) 



1 



An 1 



1 



. e -i(k t+kx) ^ dk (5 J4) 



(k + 10) - |3 2 k 2 

— oo — oo 

After the evaluation of the integral in the variable ko, G Te t takes the form: 

oo 

G ret (t,x) = -^H(t) [ e - i(p2k2t+kx) dk (5.15) 

— oo 

With a square's completion ( 15 . 1 2 j) transforms into: 



G ret (t,x) 



iHft] 



g40^t 



e ls ds 



From the result of ref . [9] 



e xsl ds = \^n:e ll < 



we have 



G ret (t,x) = -iH(t) ( 



m \ i 



Taking into account that for v = 1 : 



\2nihts 



h 



e 2m 



(5.16) 



(5.17) 



(5.18) 
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we obtain the usual form of G Tet (see ref.fTU]) 



1 



\ 2 im(x— x ) 2 

atUKt-f) ) e ^ 7T (5 ' 19) 



With a similar calculus we have for G a 



dv 



G (:dv (t-t',x-x ) = iH(t'-t) ( 2nih ™_ t] ) e (5.20) 



1 ' 2 

2 im(x— x ) 



and for Gw: 

, , i , / rn \ 2 tm(x-x ) 2 

Gw (t-t,x-x) = - 5 Sgn(t-t)( 55 j is - 7i ) (5.21) 

6 Discussion 

We have defined the fractionary Schrodinger equation for all values of the 
complex variable v and treated the cases of the free particle and the po- 
tential well. We have performed this task based in a earlier paper (ref. [5] 
where we have shown the existence of a general fractional calculus defined 
via tempered ultradistributions. In fact, all ultradistributions provide in- 
tegrands that are analytic functions along the integration path. Therefore 
these properties show that tempered ultradistributions provide an appropri- 
ate framework for applications to fractional calculus. With the use of this 
calculus we have generalized in the present work the results obtained by 
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Naber (ref.[3]). As special examples, for v = 1 the results obtained coincide 
with the usual Quantum Mechanics, and the cases v = 1/2 and v = 2 have 
shown the appearance of extra terms, besides to those with the usual (*v = 1 ) 
framework. In addition we have obtained a general expression for the Green 
function of the free particle, and shown that for v = 1 this Green function 
coincide with the obtained in ref . [TU] . For the benefit of the reader we give 
in this paper two Appendixes with the main characteristics of n-dimensional 
tempered ultradistributions and their Fourier anti-transformed distributions 
of the exponential type, and a third Appendix about the general fractional 
calculus defined via the use of tempered ultradistributions. 
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7 Appendix I: Distributions of Exponential 
Type 

For the sake of the reader we shall present a brief description of the principal 
properties of Tempered Ultradistributions. 

Notations. The notations are almost textually taken from ref[T2]. Let 
M n (res. C n ) be the real (resp. complex) n-dimensional space whose points 
are denoted by x = (xi,x 2 , ...,x n ) (resp z = (zi,Z2, ...,Zn.)). We shall use the 
notations: 

(i) x + y = (xi +t(i,X2 +\)2, — ,Xn + -y u ) ; ax = (ax b ax 2 , ax n ) 

(ii) x ^ means Xi ^ 0,X2 ^ 0, ...,x n ^ 

n 

(iii) x -y = Y_ Xjijj 

j=i 

(iV)|x|=f |xj | 

Let N n be the set of n-tuples of natural numbers. If p G N n , then p = 
(pi,P2, ...,p n ), and pj is a natural number, 1 ^ j ^ n. p + q denote (pi + 
qi>P2 + q2,...,Pn+ qn) and p ^ q means p! ^ q b p 2 ^ q 2 , ...,Pn ^ q n - * p 

n 

means x Pl x P2 ...x Pn . We shall denote by | p |= Y. Pj an d by D p we denote the 
differential operator 3 p,+P2+ - +Pn /3x 1 Pl 3x 2 P2 ...3x n Pn 

For any natural k we define x k = x k x k ...x k and 3 k /3x k = 3 nk /3x k 3x k ...3x k 



23 



The space "H of test functions such that e p|x '|D q (()(x)| is bounded for any 
p and q is defined ( ref.[12j ) by means of the countably set of norms: 

||$|| = sup e? M |D q $(x)| , p = 0,1,2,... (7.1) 

0<q<p,x 

According to reference [14"] H. is a 7C{M. p } space with: 

M p (x) = e ( P- 1)M , p = 1,2,... (7.2) 

/C{e^ p J' x '} satisfies condition (AO of Guelfand ( ref. [13J ). It is a countable 
Hilbert and nuclear space: 

oo 

K{e (v - rM } = U = f)U v (7.3) 
P =i 

where H p is obtained by completing H with the norm induced by the scalar 
product: 



v 



,2(p-1)|x| 



P 

q=0 



^D q $(x)D q {p(x) dx ; p = 1,2,... (7.4) 



where dx = dx] dx 2 ...dx n 

If we take the usual scalar product: 



OO 



4>(x)iKx) dx (7.5) 



then "W, completed with (17. 5)) , is the Hilbert space H of square integrable 
functions. 
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The space of continuous linear functionals denned on "H is the space Aoo 
of the distributions of the exponential type ( ref.[12] ). 

The "nested space" 

H=(^,H,A 00 ) (7.6) 

is a Guelfand's triplet ( or a Rigged Hilbert space [T3] )■ 

In addition we have: 'Hc«ScHc«S'c Aoo, where S is the Schwartz 
space of rapidly decreasing test functions (ref [15] ) . 

Any Guelfand's triplet G = (<D,H, O ) has the fundamental property 
that a linear and symmetric operator on <D, admitting an extension to a 
self-adjoint operator in H, has a complete set of generalized eigen-functions 
in O' with real eigenvalues. 



8 Appendix II: Tempered Ultradistributions 

The Fourier transform of a function (() G H. is 

$(x) e iz ' x dx (8.1) 



oo 
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4>(z) is entire analytic and rapidly decreasing on straight lines parallel to the 
real axis. We shall call Sj the set of all such functions. 

fj=F{H} (8.2) 
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It is a S{M. p } space ( ref. [H] ), countably normed and complete, with: 

M v (z) = ()+\z\y (8.3) 
9) is also a nuclear space with norms: 

||cb|| =sup(1+|z|) p |4>(z)| (8.4) 

zev n 

where V k = {z = (zi, Z2, Zn) G C n :| Imzj |^ k, 1 = n ) 
We can define the usual scalar product: 



< 4>(z),iMz) >= 



(j>(z)\|Ji(z) dz = 



CO 



4)(x)\Jj(x) dx (8.5) 



where: 



t|)(x) e- izx dx 



and dz = dzi dz2...dZn. 

By completing with the norm induced by (18. 5p we get the Hilbert space 
of square integrable functions. 

The dual of fj is the space U of tempered ultradistributions ( ref. [12] ). 
In other words, a tempered ultradistribution is a continuous linear functional 
defined on the space $) of entire functions rapidly decreasing on straight lines 
parallel to the real axis. 

The set U = ($), H, U) is also a Guelfand's triplet. 
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Moreover, we have: i]C5cHc5 CU. 

Li can also be characterized in the following way ( ref.[12] ): let be 
the space of all functions F(z) such that: 

I- F(z) is analytic for {z G C n : |Im(zi)| > p, |Im(z 2 )| > p, |Im(z n )| > 

Pi- 
ll Ff z)/z p is bounded continuous in {z G C n : |Im(zi)| > p, |Im(z 2 )| ^ 

p, |lm(zn.)| ^ p}, where p = 0, 1 , 2, ... depends on F(z). 

Let FT be the set of all z-dependent pseudo-polynomials, z G C n . Then 

U is the quotient space: 

III u = Ao/n 

By a pseudo-polynomial we understand a function of z of the form 

^ s Z-G(zi,...,Zj_i,Zj+i,...,Zn.) with G(zi,...,Zj_i,Zj+i,...,Zn.) G Ad 

Due to these properties it is possible to represent any ultradistribution 
as ( ref.[T2] ): 

F((j>) =< F(z), cj)(z) >= oF(z)cj)(z) dz (8.6) 

r 

r = Ti U Fz U ...r n where the path F, runs parallel to the real axis from — oo 
to oo for Im(zj) > C, C > P and back from oo to — oo for Im(Zj) < — C, 
— C < —p. ( F surrounds all the singularities of F(z) ). 

Formula (18. 6p will be our fundamental representation for a tempered ul- 
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tradistribution. Sometimes use will be made of "Dirac formula" for ultradis- 



tributions ( ref. |TT] ): 



Ffz] 



(2m) 



f(t) 



(ti -Zi)(t 2 -Z 2 )...(t n -Zn 



dt 



•7) 



where the "density" f (t) is such that 



()F(z)cj)(z) dz 



f(t)(J)(t) dt 



(8.8) 



While F(z) is analytic on V, the density f(t) is in general singular, so that 
the r.h.s. of (18. 8 1) should be interpreted in the sense of distribution theory. 

Another important property of the analytic representation is the fact that 
on P, F(z) is bounded by a power of z ( ref.jl2] ): 



|F(z)| < C|z| p 



19) 



where C and p depend on F. 

The representation (18 .6p implies that the addition of a pseudo-polynomial 
P(z) to F(z) do not alter the ultradistribution: 

o{F(z) +P(z)}4>(z) dz = oF(z)cMz) dz + oP(z)4>(z) dz 



But: 



oP(z)c|)(z) dz = 
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as P(z)4>(z) is entire analytic in some of the variables z, ( and rapidly de- 



creasing 



(){F(z) + P(z)}4>(z) dz = oF(z)cj)(z) dz 



(8.10) 



9 Appendix III: Fractional Calculus 



The purpose of this sections is to introduce definition of fractional derivation 



and integration given in ref. [11] . This definition unifies the notion of integral 
and derivative in one only operation. Let f (x) a distribution of exponential 
type and F(fl) the complex Fourier transformed Tempered Ultradistribution. 
Then: 



F(k) = H[3(k)] 



f(x)e ikx dx-H[-3(k)] 



f(x)e ikx dx 



(9.1) 



(H(x) is the Heaviside step function) and 



1 



f( x ) = _(bF(k)e- ikx dk 
in 



(9.2) 



where the contour F surround all singularities of F(k) and runs parallel to 
real axis from — oo to oo above the real axis and from oo to — oo below the 
real axis. According to [11] the fractional derivative of f (x) is given by 



d A f(x) 
dx A 



1 

2rt 



d)(-ik) A F(k)e- ikx dk + cb(-ik) A a(k)e- lfcx dk (9.3) 



—ikx 
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Where a(k) is entire analytic and rapidly decreasing. If A = — 1, d A /dx A is 
the inverse of the derivative (an integration). In this case the second term of 
the right side of (19. 3 p gives a primitive of f(x). Using Cauchy's theorem the 
additional term is 

l^-e - ikx dk = 27ra(0) (9.4) 
k 

Of course, an integration should give a primitive plus an arbitrary constant. 
Analogously when A = —2 (a double iterated integration) we have 

o-^e- ikx dk = y + 5x (9.5) 

where y and 6 are arbitrary constants. 
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